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Abstract. We study the Casimir effect for scalar fields with general curvature coupling subject 
to mixed boundary conditions (1 + /3 m ra M <9 M )(/9 = at x = a m on one (m = 1) and two (m = 1, 2) 
parallel plates at a distance a = a<i — a\ from each other. Making use of the generalized Abel- 
Plana formula previously established by one of the authors ]l]], the Casimir energy densities 
are obtained as functions of /3i and of /3i,/?2,a, respectively. In the case of two parallel plates, a 
\ decomposition of the total Casimir energy into volumic and superficial contributions is provided. 

The possibility of finding a vanishing energy for particular parameter choices is shown, and the 
existence of a minimum to the surface part is also observed. We show that there is a region 
in the space of parameters defining the boundary conditions in which the Casimir forces are 
repulsive for small distances and attractive for large distances. This yields to an interesting 
■ possibility for stabilizing the distance between the plates by using the vacuum forces. 

o 

o . 

1 Introduction 

Although the existing literature about the Casimir effect is quite sizable in volume (for reviews 
see. e.g. [§]), we feel that relatively little attention has been devoted to quantum fields subject 
to Robin — or mixed — boundary conditions on plates. A possible reason is that this type of 
condition appears when decomposing the modes of the electromagnetic field in the presence of 
perfectly conducting spheres (see refs. |3||-[f§), but are not required in the analogous problem 
with parallel plates, where the mode set can be divided into eigenmodes satisfying Dirichlet and 
Neumann conditions separately. 

However, Robin conditions can be made conformally invariant, while purely-Neumann ones 
cannot. Thus, Robin-type conditions are needed when one deals with conformally invariant 
theories in the presence of boundaries and wishes to preserve this invariance. The importance 
of conformal invariance in problems related to the Casimir effect has been emphasized, e.g. in 
refs. |], (see also |^). The Casimir energy- momentum tensor on background of conformally 
flat geometries can be obtained by the standard transformation from the corresponding flat 
spacetime result (see ||). For instance, by this way in [|l(J the Casimir effect for a scalar 
field with Dirichlet boundary condition is investigated on background of a static domain wall 
geometry. To derive the Casimir characteristics via conformal transformations for the case of 
Neumann boundary condition we need to have the corresponding flat spacetime counterpart 
with Robin boundary conditions and Robin coefficient related to the conformal factor. It is 
interesting to note that the quantum scalar field satisfying Robin condition on the boundary 
of cavity violates the Bekenstein's entropy-to-energy bound near certain points in the space of 



X 

5-H 



the parameter defining the boundary condition [11]. The Robin boundary conditions are an 



extension of the ones imposed on perfectly conducting boundaries and may, in some geometries, 
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be useful for depicting the finite penetration of the field into the boundary with the "skin-depth" 
parameter related to the Robin coefficient | |12| |. On the other hand, the relevance of mixed-type 
boundary conditions to spacetime models and quantum gravity has been highlighted in refs. 
[pHy , [p~I|], This type of conditions naturally arises for the scalar and fermion bulk fields in the 
Randall-Sundrum model fl5|| . In this model the bulk geometry is a slice of anti-de Sitter space 
and the corresponding Robin coefficient is related to the curvature scale of this space. 

In the present work we discuss several aspects of the Casimir effect for a massless scalar 
field, with curvature coupling, obeying Robin boundary conditions on one or two parallel plates 
on background of a flat P-dimensional spacetime. The dimensional dependence of physical 
quantities in the Casimir effect is of considerable interest and is investigated for various types 
of geometries (see, for instance, 1C], [17], ]l8|). In sec. § we explain how Robin conditions can 
adopt a conformally invariant form. The Casimir effect with one plane boundary is considered 
in sec. [3|, while sec. [| is dedicated to the set-up where two parallel plates are present. Then, the 
volume and surface contributions to the total Casimir energy (for this second case) are analyzed 
in sec. ||. Our ending comments follow in sec. |6|. 



2 Conformal invariance and boundary conditions 

Let's consider a massless scalar field if with curvature coupling £ on background of a T>- 
dimensional spacetime manifold M with boundary dM. The action for this field is 



d v x 



-9<pp + ZR\ip, 



(2.1) 



M 



where □ - is the Laplace-Beltrami operator. The Lagrangian corresponding to ( |2.1[ ) differs 
from the often used Lagrangian by a total divergence leading to the additional surface term 
\ fdM d Dlx V~9 n^tpd/if with being the unit normal vector to dM. As it has been noted in 
ref.[0], this term plays a crucial role in the cancellations between surface and volume divergences. 
Note that the additional surface term is zero for Dirichlet and Neumann boundary conditions 
on dM, but is nonzero for the more general Robin case. 

Consider a conformal transformation realized by a Weyl rescaling of the spacetime metric 



9fiu(x) — ► VL 2 {x)g lJiV (x). 
Under these transformations, the (p field will change by a rule of the type 

(p(x) — > Q a (x)tp(x). 
As a result, the action undergoes the following transformation: 



(2.2) 
(2.3) 



s[n a ip,n 2 g] 



1 



X>-2+2a 



□o 



<pp + CR]<p + [a + 2£(P - 1)]-^V 



+ (2a + V-2)<r%£.<pd v <p 
+[a(a + T>-3)-£(p- 1)(V - 4)] g 



(2.4) 



The action S will be invariant if T> — 2 + 2a = and all the terms containing derivatives of 
vanish. These two requirements are satisfied provided that 



V 



a 



2 1 
V-2 

4(2? - 1) 



(2.5) 
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Next, we shall consider the effect of the transformation on a boundary condition of the Neumann 
type 

n! i V ll <p(x) = 0, (2.6) 

where is a normal space-like vector (i.e., g iiv n l> "n v = — 1) perpendicular to the boundary, 
and covariant derivative reduces, in this case, to the ordinary partial derivative because ip 
is just a scalar function. Let denotes the transformed version of n^. If we require that the 
normalization be maintained, we shall have £l 2 g ^n^rf = — 1, whose solution is 

n» = ^n". (2.7) 



Taking into account ( p. 2D , (|2.3|), ( |2.5| ) and ( |2.7| ), we realize that the l.h.s of the boundary 
condition ([2.6j) transforms as 



n"V„p — » n~— InPVpip . (2.8) 

The presence of the second term indicates that a boundary condition of purely-Neumann type 
cannot be maintained under general conformal transformations. 

Similarly, if, instead of (|2.6|), one takes a generic Robin boundary condition 



n"V„)p(z)=0, (2.9) 
one can readily observe that it changes according to the rule 

(^ + n"V M )^^Jr— InVcp + n^^ ^-<p\ , (2.10) 

where ^ indicates the result of transforming the function. The boundary condition ( |2 .9[) 

can be preserved only if the transformed version is proportional to the initial form. Thus, one 

_ — . p— i 

demands that the r.h.s of ( [2.10D be equal to Q ~ (^ + n^V^) (p. This leads to a specific 

transformation rule for the Vl/ function, which reads 

as already observed in ref.0. 

Now, suppose that we have a valid ^ function satisfying ( 2.11| ). We can consider n along 
the x-axis and set boundary conditions on the planes x = a\ and x = 02. Provided that 
^>(x = ai) 7^ and *$>(x = 02) 7^ 0, we may write the boundary conditions at these points in the 
form 

{l + p m {-l) m - 1 d x )p = atx = a m , m = l,2, (2.12) 

where 

P™ = *n-^ — 7, ™ = 1 > 2 - (2.13) 

V(x = a m ) 

One may consider the subgroup of transformations in which does not depend on the x- 
coordinate. Then, a possible ^ is given by *[(<?)] = (—g) 1 ^ 2T '- These particular transformations 
correspond to the restriction of the initial group to planes parallel to the plates. In a strictly 
Euclidean or Minkowskian spacetime, this form of Vl/ would imply 02 = f3\. 
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3 Casimir stresses for a single plate geometry 

In this section we will consider scalar field in D-spatial dimensions — thus, D = D — 1 — with 
general coupling £ satisfying Robin boundary condition on the single boundary x = on back- 



ground of a flat spacetime. Such a situation is like limiting eq. (|2.12|) to m = 1 only, and with 
a\ = 0, i.e., 

(l + ftn"5Xt,x) = (l + ft^Mt.x) =0, x = 0. (3.1) 

Here we consider the vacuum fluctuations in the region x > 0. For the region x < the boundary 
condition has the form (1 — /?i d x )<p(t, x) = at x = 0. The corresponding results can be obtained 
from the previous case replacing f3\ — ► — 

The eigenfunctions satisfying boundary condition ( |3.1| ) are in form 

<^ k (t, x) = ====== cos(£;x + ai), (3.2) 

V 2- u ~ i 7r- u u; 

where k = (fc, kj_), a; = |k|, < k < oo, and 

smoi = — =====, cosai = — , 3.3 

In the case 0% > there is also a purely imaginary eigenvalue = with the normalized 
eigenfunction 

where k± = |kj_| > 

3.1 Vacuum densities 

From the symmetry of the problem it follows that the vacuum expectation values (v.e.v.) for 
the energy-momentum tensor (EMT) have the form 

<0|2*|0) = diag(e, -p, -p ± ,..., - P± ). (3.5) 

The corresponding energy density e and effective pressures p, p± can be derived by evaluating 
the mode sum 

(0\T ik (x)\0) = Y,TikWa(x),^ a (x)}, (3-6) 

a 

where (p a (%) = (vh.^Vk^j i s a complete set of solutions to the field equation and the bilinear 
form Tik{f,g} is determined by classical energy-momentum tensor for the scalar field (see, e.g., 
||). Using the field equation we will present it in the form 



T ik = di<pd k (p + 



^ (3.7) 



Using formula ( |3.6| ) with eigenfunctions ( |3.2| ) , ( |3.4| ) and EMT from ( |3.7| ) one finds (no summation 

over i) 



(2?r)D J dk^ 1 cos (2A;x + 2ai) + 
2ft 7 (2vr)^ Jk? ± -1/P* 
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where 9(x) stands for Heaviside step function, 
A (k) = k 2 ± + iCk 2 , A 1 (k) = 0, Mk) 



D - 1 



(4f - l)k 2 



2,...,D, 



and 



(0 M \T lk (x)\0 M ) 



(2ir) D 



00 dk 
^ 



diag ( u> 2 , k 2 , 



k\ 



D-V"' D-l 



(3.9) 



(3.10) 



are the corresponding quantities for the Minkowski vacuum \0m)- The last summand on the 
right of ( |3.8[ ) corresponds to the contribution of eigenfunctions ( |3.4[ ) . In this term the integration 
over kj_ goes in the region |kj_| > l/d\. 

In eq. (3.8) the integrals over kj_ may be evaluated by using the formulae 



d D - l k L g(k ± 

x y n dy 



2n (D-l)/2 j-oo 

W) 



1 



vV + 1 

00 y n dy 



Vy 2 ^- 



20F 
l 



k?- 2 g(k ± )dk A 
n + V 



1 / 71 + 1 \ nn 



(3.11) 
(3.12) 
(3.13) 



As a result for the difference between v.e.v. for vacuums |0) and |0m) one obtains 
(TiklsUB 



(0|T ifc (x)|0) - (0 M |T ife (x)|0 M ) = {ic ~ 2 ^^g/ 2 + 1 diag(l, 0, -1, . 



dkk D cos (2kx + 2a>i{k)) + 2tt0(/3 1 



Hi 



ttD 

d+t cos — 



-i) 

(3.14) 



where ^ = £ c = (D — 1)/AD corresponds to the conformally coupled scalar field (recall eq. ( |2.5[) 
with T> = D + 1) and the function ot\{k) is defined as (3.3). Using these relations the integral 
over k can be presented in the form 



Jo 



dkk D cos (2kx + 2a x {k)) 



dkk D cos (2kx) 



00 „ , n cos (2kx) + kBi sm(2kx) 
dkk 



k 2 3\ + 1 



(3.15) 



To evaluate the second integral on the right we will use the formula 



00 k n ~ 1 e~ p ° k 
k 2 + z 2 



v n-2 



dk 



r(n - 1) r e w+™/2 r ( 2 



z) + e -*»*-«Trf/a r (2 _ U) _ ipoZ ) 



(3.16) 



where Re (z), Re (n), Re (po) > 0, and r(a,z) is the incomplete gamma-function (see, e.g., [|20|). 
The corresponding integrals in ( 3.15| ) with cos and sin can be obtained as real and imaginary 
parts with po = a — 2ix, a —> +0. Using the formulae [pcf ] 



T(-n,y) 



(-iy 



in 



n-l 



(-l) J j! 



3=0 



E 1 (-y±i0) = -Ei(y)T«7r 



(3.17) 



(3.18) 



the corresponding v.e.v. can be presented in the terms of the functions E\(y) and Ei(y). As a 
result it can be seen that the difference between v.e.v. on the left of ( |3~l4T ) is finite for x > 
and can be presented in the form 
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00/ SUB 



4D(e-6) 



7r^/ 2 r( J D/2 + l)(4|x|) D - 



D+l 



-IvlEJif 



D-l 



and 



00/5(75 ' 



, y = 2x^1 



(3.19) 



(3.20) 



(Tu) SUB — 
( T u)sUB = °- 

Note that in these formulae the term in the v.e.v. (|3.14| ) coming from the contribution of the 
eigenfunctions ( |3.4| ) (and divergent for even values of D) is cancelled with the corresponding 
term coming from the integral ( |3.15| ). As a result the subtracted v.e.v. are finite for all values 
y>0. 

We have considered the region x > 0. As it have been mentioned above the corresponding 
formulae for x < can be obtained from ( |3.19| ) replacing j3\ — ► — j3i. As a result the value of 
y remains the same and the distribution for v.e.v. ( 3,19| ) is symmetric for regions x > and 
x < . The cases for the Dirichlet and Neumann boundary conditions are obtained from (3.19) 
in limits /3i — ► and @\ — ► oo, respectively: 

(( r 00)5[/ B )Dirichlet = -((^bo) SUb) Neumann = ^^p^j/^j D+l F ( 2 ) ' ( 3 " 21 ) 

Note that the first term in ( |3.19| ) coincides with the energy density corresponding to Neumann 
case. As follows from ( |3.19| ), the regularized v.e.v. is zero for a conformally coupled field. This 
result can be obtained also without explicit calculations by using the continuity equation and 
zero trace condition for the EMT. 

Using the asymptotic formulae for the functions E\{z) and Ei(z) from ( |3.19 ) one obtains the 
following asymptotic expansion of the vacuum energy density for large y: 



00/ SUB 



00 1 SUB) Dirichlet 



+ 



D(£ - He 



CO 



Ji 



2 2D-l 7r (D+l)/2 



I D+l yj-D- 

j=D+l y 



(3.22) 



As we see at distances far from the plate the vacuum energy density coincides with that for the 
Dirichlet case and is positive for £ > £ c and is negative for £ < £ c . At distances near the plate , 
|y| <C 1 the vacuum energy density is dominated by the first summand in the figure brackets in 
( |3.19| ), As it have been noted this summand coincides with the energy density for the Neumann 
case and hence has opposite sign compared to the case of far distances. As a result the vacuum 
energy density has a positive maximum or negative minimum (depending on the sign of £ — £ c ) 
for some intermediate value of x. Near the plate from ( 3.19j ) we have the following expansion 
for the vacuum energy density 



00/ SUB 



(( T 00>S[/B) N eumaiin 




D-j 



y D+1 In \y\ 



+ 0(1). (3.23) 



All terms presented lead to the diverging contribution to the energy density at the plate surface. 
These surface divergences are well known in quantum field theory with boundaries and are 
investigated near arbitrary shaped smooth boundary for the minimal and conformal scalar and 
electromagnetic fields H,@. 
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Figure 1: Regularized Casimir energy density for a single plate geometry — formula (|3.19 ) 
multiplied by — £c) as a function of y for D = 1 (left) and for D = 3 (right). 





yi 


2/2 


( T oo) SUB (vi) 


( T oo) sub (v-z) 


1 


2.83 


-2.54 


0.136 


0.021 


3 


4.83 


-4.53 


2.04 10" 3 


1.68 10~ 4 


5 


6.83 


-6.52 


1.19 10~ 5 


6.68 10~ 7 


7 


8.83 


-8.51 


4.44 10- 8 


1.89 10~ 9 



Table 1: Maxima and corresponding values for the energy density (multiplied by \Pi\ D+1 /(£,— £c)) 
in the case of a single plate. 

In Fig. [j] we have plotted the dependence of \(3i\ D+1 e/(£ — £ c ) on y for the cases D = 1 
and D = 3. As we see, for a given f3\ the vacuum energy density has a maximum for x = x^ = 
Piyi(D), i = 1,2, where y = y% correspond to the maxima in the figures and depend only on 
the space dimension D. In Table |] we show the values for these maxima and the corresponding 
values of the energy density for D = 1,3,5,7. 



When j3\ — > 0, one has xi — > and (Tqq) sub (x = Xi) ~ \(3\\ 
distribution corresponding to the Dirichlet boundary condition. 



+oo and we obtain the 



3.2 Total vacuum energy 

Integrating ( 3.19| ) in x, we find for the energy per unit surface area from x = x\ to x = oo, 



E Xl < 



x<oo 



dx (Too) 



SUB 



X'l 



-on - £ c 



2 d tt d / 2 F(D/2 + l)(3[ 



D 



r(g) 
y? 



D-l 



2e-»Ei(yi) + 2 J] 



(i-i)! 



3=1 2/1 



yi 



(3.24) 
2xi 



To investigate the total vacuum energy it is convenient to present the subtracted vacuum EMT 
( |3.14 ) in another alternative integral form. In the integral term on the right of eq.( 3.14 ) we 
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write the subintegrand in terms of exponential functions and rotate the integration contour by 
angle ir/2 for the term with e* fc ' x ' and by angle —tt/2 for the term with e~ ik ^ x K This procedure 
leads to the relation 



-27r6>(/?i) d x cos • 



/ dkk D cos (2kx + 2 ai (k)) = -sin — p.v. / dt t D e~ 2tlxl - (3.25) 
Jo 2 Jo Pit - 1 

^rr- cos — . 

where the second summand on the right comes from the poles ±i/(3i in the case j3\ > and 
the symbol p.v. means the principal value of the integral. Now comparing to ( 3.14 ) we see that 
this summand cancels the contribution from the mode fl3.4|). Substituting ( |3.25|) into ( p. 14 ) and 
using the gamma function reflection formula 

irD ( D\ vr , 

sin r = ; — ; (3.26) 

2 V 2 J r(D/2+l) V ; 



for the subtracted v.e.v. one finds 



{Tik) 



SUB 



2 D - l 7T D / 2 T(D/2 + 



-p.v. r^^ e -^N^±l dia g ( i,o,-l,...,-l). (3.27) 
1) Jo pit — I 



Integrating the corresponding energy density for the total volume energy in the region < x < oo 
we receive 

E< " ,M <« ^r (D /2 + i) p - v - f s£r (3 ' 28) 

Here and below the index (s) is for quantities describing the single plate geometry. 

In addition to volume energy fl3.28|) there is also surface energy contribution to the total 
vacuum energy. The corresponding energy density is defined by the relation (see, Q) 



T( gurf) = _jL-±S(o-:<nf) r y) rY ~. (3.29) 



where 5(x;dM) = 5(x + 0) is a "one sided" 5 - function. From this formula it follows that 
the surface term is zero for Dirichlet or Neumann boundary condition (as the factors <p or d x <p 
would then vanish) but yields a nonvanishing contribution for Robin boundary conditions. The 
evaluation procedure for the corresponding v.e.v. is similar to that given above for the volume 
part and leads to the expression 

(OIT^IO) = + (3-30) 



J clkk ' 1 sin (2kx + 2ai(fc)) - 2-^6 (Pi) ^ cos — 



where the second summand in the square brackets on the right comes from the purely imaginary 
mode ( |3.4p . Transforming the integral term by the way similar to (3.25) and integrating the 
energy density over the region < x < oo for the total surface energy one finds 

£ C)(.u,.) (A) _ D J-W n, v . r dff ^M±i. (3.31) 
2 D 7T D / 2 T(D/2 + 1) P J Pit -I V ' 

As in the volumic part the contribution of mode ([O]), divergent for even values of D, is cancelled 
by the term coming from the poles ±i/(3\, (3\ > 0. The total vacuum energy is the sum of the 
volumic and surface parts: 



sl *»-^pTi) p ''i" i,M ^' <3 ' 32) 
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The £ - dependence in the volumic and surface energies cancelled each other and we have £ - 
independent total vacuum energy. The expressions (3.28), ( 3.31| ) and ( |3.32| ) are divergent in 
given forms. Due to the relations 

£(s)(vol) = w ^ _ Q E ( S ) ^ E (s)(suri) = d q _ 4 ^ E (s) (3_ 33 ) 

it is sufficient to regularize the total vacuum energy. We will return to this question in section ||. 
By taking into account that the regularized value of the integral dtt 0-1 is equal to zero we 
conclude from ( |3.32| ) that the total regularized vacuum energy is zero for Dirichlet and Neumann 
scalars (see also M). 

4 Scalar Casimir effect with Robin boundary conditions on two 
parallel plates 

In this section we will consider a scalar field with £ coupling and satisfying Robin boundary 
conditions ( |2.12| ), i.e., 

(l + (3 m n»d fl Mt,x) = (l + f3 m (-l) m - 1 d x Mt,x)=0, x = a m , m = 1,2 (4.1) 

on plane boundaries x = a\ and x = a<i- The corresponding eigenfunctions in the region between 
the plates can be presented in two equivalent forms (corresponding to m = 1,2) 

<Pk(t,x) = pe ik ^- iujt cos[k\x-a m \ + a m ], (4.2) 



where k = (k, kj_), oj = < lk\ + k 2 , and 



1 k(3 m 
sma m = — =, cosa m = — =. (4.3 

Vk 2 Pl + l Vk 2 (3 2 m + 1 

From the boundary conditions one obtains that the eigenmodes for k are solutions to the fol- 
lowing equation 

F(z) = (l — 6i&2-2 2 ) sin z — (62 + b\)z cos z = 0, z = ka, bi=@i/a, a = a% — a\. (4.4) 

In Appendix [A] the values of the coefficients hi are specified for which this equation has purely 
imaginary zeros. For these zeros the region of the allowed values for k^ is restricted by the 
reality condition of uj: k\ > -k 2 . 

The coefficient (5 is determined from the orthonormality condition 

J (p k ^dV = 7^<Hk± - ^±)S k k', (4.5) 

where the integration goes over the region between the plates. Using the form of the eigenfunc- 
tions one obtains 



/? 2 wa(2vr) 



D-l 



1 + - — sin(/ca) cos(ka + 2a r 
ka 



(4.6) 



(on the class of the solutions ka = z to Q4.4| ) the expressions on the right are the same for m = 1 
and m = 2 ). 



9 



4.1 Vacuum energy density and stresses 

The v.e.v. for the EMT can be found by evaluating mode sum ( |3,6| ) with the energy-momentum 
tensor ( |3.7| ). Substituting the eigenfunctions (|4.2| ) for the vacuum EMT components one finds 



$ik sr^ f d D 1 k^ Bi + Ai(z/a) cos [2z \x - a m \ /a + 2a m ] 



(0\T ik (x)\0) = ^ ^ / "v., — — """" ' ^ m J , (4.7) 
2 z= t^ ■/ W Jk 2 ± a 2 + z 2 [1 + \ sin 2 cos(z + 2ai)] 



where 



and the coefficients Aj(&) are defined as ( |3.9| ) with k = z/a, z = X n are positive solutions to 
equation (fO|), and z = iyi, yi > are the purely imaginary solutions in the upper half-plane. 



k 2 

B = k 2 ± + z 2 /a 2 , B 1 = z 2 /a 2 , B 2 = i = 2, (4.8) 



For the latter case the kj_-integration in (4.7) extends over the region > yf/a. 

Next, we apply to the sum over n in eq. (|4.7[) the formula derived in the appendix [b[ Note 
that ct\ depends on z = X n ,iyi as well and 

f n \x-a m \ \ z 2 b 2 m -l f x-a m \ 2zb m . / \x - a m \\ . , 

cos (2*— ^— + 2a m j = cos (2*-^- j - sm (2*— ^— j . (4.9) 

In ([4.7]), we perform the integration over by using formula ( |3.11| ). Further, introducing a new 
integration variable y = k±/k and evaluating the corresponding integrals over y using formulae 
(fn|) and flOp , we find that the vacuum EMT has the form ([3.5|). Energy density, e, pressures 
in perpendicular, p, and parallel, p±, to the plates directions are determined by relations 

q ( x > = " np+i d/2 D+i 1^ 1,1- / i o \ ? 9 = e,P,P±, (4.10) 

^ 1 + - smzcoslz + 2ai) 

z=\ n ,±i yi z 

where the prime on the summation sign means that the contribution of the terms corresponding 
to the purely imaginary zeros have to be halved and the following notations are introduced 



z 2 b 2 m + 1 



(z 2 b 2 m - 1) cos ( 2z X -^ ) - 2zb m sin ( 2z l * 



a J \ a 

f£>(z,x) = Dz D , 

fm ( z ,x) = — fm {z,x). 

(4.11) 

It follows from here that p± = — e, and, for the conformally coupled scalar (i.e. £ = £ c ), the 
components of the vacuum EMT are uniform between the plates. Similarly to the cases of 
Dirichlet and Neumann boundary conditions, these properties can be also directly obtained by 
using symmetry arguments. The field equation and boundary conditions are invariant with 
respect to the Lorentz boosts in directions parallel to the plates. It follows from here that the 
corresponding (transverse to x axis) part of the vacuum EMT is proportional to the metric 
tensor, and hence 

e = (0\T°(x)\0) = (0\T 2 (x)\0) = ... = (0|T^(x)|0) = -p ± . (4.12) 

For the conformally invariant field, from the zero trace condition one finds p = —(0\T^(x)\0) = 
De. By the symmetry of the problem, the quantities q depend only on x coordinate. From 
the continuity equation for the EMT it follows that p'{x) = 0, and therefore the vacuum EMT 
is constant. Unlike the case of Dirichlet and Neumann boundary conditions, for the confor- 
mally coupled scalar the functional dependence on the plates separation cannot be determined 



10 



by purely dimensional arguments, because we have three parameters with length dimension, 
a,(3i,^2- To obtain the dependence on these parameters we need an explicit calculation. 



The summation in ( 4.10| ) can be done by using formula 
result one obtains 



7Tfm\z,x) 



1+ shp cos ( z + 2 , 
z=\„,±iyi 2 v ' 



f$ (z, x)dz — 2 sin 



+8.cos(- 



ttD 

-e (b 



taking f(z) = fm (z,x). As a 

F^\t,x)dt 

(6it-l)(6at-l) of. _ i 
(M+1)(M+1) 



p.v. 



exp 



ab r , 



, (4.13) 



where the following notations have been introduced 



Frh ^ (-2; x) 



1 + 



(z,a?) 



-Fi £) (z,x). 



1 



(z 2 i& + 1) cosh ( 2 



x — a r 



2zb m sinh ( 2z 



F (P) 

1 m 

Fm (z, x) 

(4.14) 

Note that on the right of ( |4.13| ) we have included the term ( B.8 ) coming from the poles zizi/b m 
for b m > . 

By comparing to the formula ( |3.14j ) it can be easily seen that the contribution of the first 
and last summands on the right of formula ( [4.13| ) corresponds to the vacuum EMT components 
for the geometry of a single plate placed at x = a m . This can be also seen by taking the limit 
(— l) m a m i — > oo, m! m, m! = 1,2, when the second integral on the right of eq. (4.13), 
multiplied by a~ D_1 , tends to zero. The regularization for the single plate case was carried out 
in previous section. Hence, using formula fl4.13| ) for the regularized v.e.v. of the EMT for the 
case of two plate geometry from ( 4.10 ) one obtains 

regq(x) = iegq$(x) 



a D + l T{D/2 + l 



-p.v. 



F$(t,x)dt 



(M-l)(M-l) 2t 



1 



(4.15) 



(6 1 t+l)(6 2 t+l)' 

Here we have used the gamma function reflection formula ( ^.26 ). In ( [4.15 ) the term reg qm(x) 
is the regularized v.e.v. for the case of a single plate placed at x = a m . This geometry was 
investigated in the previous section and the corresponding regularized quantities are given by 
relations ( |3.19| ) and (3.20) with the replacement As follows from ( [4.14j ) and ( |4.15|) 

the vacuum perpendicular pressure, regp, is uniform in the region between the plates: 



reg p 



2- D ir- D / 2 
'a D + l T(D/2 + 1 



-p.v. 



t D dt 



(M-iXM-i) o2 t 



(M+lKM+l) 



1 



(4.16) 



The Casimir forces per unit area on the plates are equal to this vacuum pressure and are 
attractive for negative values of reg p, and repulsive for positive values. In the limit |/?i|/a <C 1 
using the value of the integral 



f 

Jo 



t u dt Cr(D + 1) 



,2t 



1 



2 D+l 



T{D + 1) 



{C,r meaning the Riemann zeta function) from ( 4.16 ) one has 



reg p = (reg £>) D irichiet 



1 + 



(£) + !) + ••■ 



(4.17) 



(4.18) 
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where 



(reg p) 



Dirichlct 



D( R (D + 1) 



D + l 



(4.19) 



( 47r )(i5+i)/2 a D+i V 2 

is the corresponding vacuum pressure for the case of the Dirichlet boundary condition The 
latter corresponds to the attractive vacuum force for any value of D. Similarly, in the limit 
\(3\\/a 3> 1 one derives 



reg p = (reg p) Neumann 



1 + 4a 



1 1\ Cr(d-i) 
Pi 02 J D( R {D + 1) 



+ 



(4.20) 



with (reg p)Neumann = (reg p)Dirichiet being the vacuum pressure for the Neumann scalar. For 
given Robin coefficients Pi formula ([4 . 1 8|) gives the asymptotic behaviour of the vacuum forces 



for large values a, and eq. ( 4.20 ) gives the asymptotic for small distances. In both these limits 
the vacuum forces are attractive. 

Let us consider in detail the special case Pi = (Dirichlet condition on the plate x = a\). 



Now in the limit of large distances, \p 2 \/a 1, the vacuum forces are dominated by term ( 4.19 ) 
and, hence, are attractive. At small distances, |/?2|/o 1, from ( f4,16| ) with fix = one has 



reg p 



(reg p) 



Dirichlet 



[l + 0(a/\p 2 \)) 



(4.21) 



which corresponds to the repulsive force. It follows from here that at some intermediate value of 
a the Casimir force vanishes. As a result we have an example when the vacuum forces stabilize 
the distance between the plates. This feature is seen from Fig.||, where the Casimir force is 
plotted versus ajP% in the case D = 3 and Pi = 0. 



CD 

u 



0.5 



-0.5 - 




Figure 2: Casimir force per unit surface multiplied by 10 4 /3|, i.e. 10 4 /5|reg p as a function of 
b~ l = a/fo, for D = 3, ft = 



In ( [4.15 ) taking the limit (— l) m a m — ► oo one obtains integral form ( [3.27 ) for the v.e.v. in 
the case of a single plate. By using this formula the corresponding quantities ( 4.15| ) for two 
plate geometry can be presented as 



vegq(x) = vegq( s \x; a±) + iegq^ s \x; a 2 ) + Aq(x), 
where "interference " term has the form 



Aq(x) 



a^r( J D/2 + l) P - V -7o fe^jj^ eB-i 

(M+l)(M+l) 



AF^(t,x)dt 



(4.22) 



(4.23) 
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with notation 



AF^ e \t,x)=t D 



1 + 2D{H - £ c ) £ ^4^4 exp ( -2t lX " ""' ! 



m=l 



b m t - 1 



(4.24) 



In ( |4.22|) the quantities Aq(x) are finite for all values a\ < x < a 2 and the surface divergences 
at contained in the summands vegq^(x;a m ). Note that the quantities ( f4.23| ) are 

symmetric under the replacement b\ b 2 - 

For the conformally coupled scalar field reg<7^(a;) = 0, and for the vacuum EMT components 
between the plates from ( 4.15Q one obtains 



reg e 
regp ± 



reg p/D = e 
-reg e, 



(i) 

c 



(4.25) 



where is determined from ( |4.16| ). Hence is the vacuum energy density between the plates 
for the conformally coupled scalar field. 

In the case b\ = —b 2 using the value of integral (4.17) and the duplication formula 



T(2z) 



2 2 *- 1 ( 1 



eq.( 4.15| ) turns into 



rege 



Oe(/^+l) r fl ! 



(4.26) 



(4.27) 



^)(D+i)/2 a D+i- y 2 

This formula coincides with that derived in []2l| for the Dirichlet and Neumann cases. Notice 
that, whenever b\ = —b 2 = b, the particular value of b does not matter, as one might have been 
observed at the beginning, from the form of the F function. 



4.2 Total volume energy 



From relation (4.22) for the total volume energy per unit surface one has 



ra,2 

£( vo1 ) = / rege = £ (s) (ai < x < a 2 ;/?i) + £ (s) (ai < x < a 2 -(3 2 ) + AE, 

J ai 



(4.28) 



where E^ (a\ < x < a 2 ; m ) is the vacuum energy in the region a\ < x < 02 due to a single plate 
at x = a m with Robin boundary condition (|4.1| ). The "interference" term AE = dx Ae(x) 
can be presented in the form 



AE 



2- D TT- D / 2 

a D T(D/2 + 1 



-p.v. 



t D dt 



(bit- l)(ft 2 /.-I )„<;/. 
(M+1)(M+1) 



1 + m - tc 



e ~ 2t A b m t + 1 



E 

m=l 



bmt - 1 



(4.29) 

The energy < x < a 2 ; 0i ) can be obtained subtracting from the total volume energy of 

a single plate (given by ( 3.28[ )) the vacuum energy in the region a 2 < x < 00. The latter can be 
easily derived integrating the 00-component of eq. ( 3.27 ). As a result one finds 



E^iai <x<a 2 ;(3 1 ) = E^^Xfa) + 



2 d tt d / 2 T(D/2 + l)a D 



p.v. 



~dtt D ^e~ 2t(3lt + ] 



(4.30) 
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Taking into account this and similar formula for £A s )(ai < x < 02; (3 2 ) from flOg) and (|4~29|) 
we receive 



£(v°D = ^ ^( s )( vol )(/3 m ) + a (e« + 4D(£ - &)e. 



'C / ) 



(4.31) 



m=l 



where defined as in ( |4.16| ), and we have introduced notation 



,(2) 



61+62 



2 D TT D l 2 a D + 1 T 



t D (l - hb 2 t 2 



prr- D.V. / dt 

l + # Jo 



L t) 2 (l - M) 2 e 2 * - (1 - 6f* 2 )(l - 6|t 2 



-. (4.32) 



For Dirichlet (b\ = 62 = 0) and Neumann (b\ = 62 = 00) scalars this term vanishes. Note that 
the volume energy ( |4.31| ) is symmetric under the replacement b\ ^62- 



5 Total Casimir energy 

Up to now, we have obtained local energy densities from expectation values of the energy- 
momentum tensor. Here, we will take the alternative approach of calculating the integrated 
Casimir energy per unit volume — e c — from the eigenvalue sum, namely, 



2a 



(5-1) 



where, k n = z n /a, being {z n = \ n ,iy n } the set of the zeros of the F(z) function defined by 
eq.(|4.4|). In the case of the purely imaginary zeros the k^-integration ranges over the region 



k\ > y 2 Ja 2 . 



5.1 Zeta function regularization 



As it stands, the r.h.s. of ( |5.1| ) clearly diverges, but we shall evaluate by zeta function regular- 
ization, i.e., turning ( |5,lD into the function 



£ c (/i;s) 



J± 

2a 



d 



D-U 



(2*) 



D-l 



k 1 ~\~ k. 



2 \ -8/2 

n 



(5.2) 



and adopting the prescription that the regularized value of e c will be £ c (//; s = —1) (some of the 
pionnering works in this sort of technique are listed as ref. p2}| ). Note that we have introduced 
the arbitrary mass scale \i in order to keep dimensionless the quantity raised to the power of 
—s/2. After integrating over kj_, eq. (|5.2|) may be written as follows 



e c (/i;s) 



1+1 



r(f) 



2 D 7T (D-i)/2 a D - s T(l) 



CA(a) 



cos(vrD/2) 
sin(7rs/2) 



where 



a = s + l-D, (5.3) 



(5.4) 



71=1 



is the zeta function for the real zeros of F(z). 

First, we look at the asymptotic form of the A n 's in order to find the convergence boundary of 
Ca((t). If 6162 / 0, then F(z) ~ b^z 2 svnz as \z\ » 1, and Ca(c) ~ ^(7rn) _<T = ^"Cr^) — Cr 

n>l 

denoting the Riemann zeta function — , which has a pole at a = 1. When b\ = and 62 7^ (or 
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&i ^ and b 2 = 0), one has F(z) ~ — 6 2 ,i-zcos2: for \z\ > 1. Thus, Ca(ct) ~ ^(vr(n + l/2))- <T = 

7r~ cr (2 CT — 1)£r(ct) which has a pole at a = 1, too. Finally, if &i = 62 = 0, = sinz and, 

therefore, Ca( (7 ) = 7r_ °'G?( cr ) exactly. In view of this, we realize that Ca(o") has its rightmost 
pole at a = 1, regardless of the values of 61, b 2 . 

This way, we see that eq. (5^) is, initially, just valid for the domain Re a > 1. Nevertheless, 
in order to obtain the regularized Casimir energy, one has to find the analytic continuation of 
( |5.4D to a = —D, which corresponds to s = — 1 (here, D = 1, 2, 3, . . .). This task will be done 
by analytic extension of an adequate contour integration in the complex plane. An immediate 
consequence of the Cauchy formula for the residues of a complex function is the expression 



Ca(ct) 



— / dzz- a — lnF(z), forReo->l, 



2m 



c 



dz 



(5.5) 



where C is a closed circuit enclosing all the zeros of F(z). In this case, we assume that C is 
made of a large semicircle — with radius tending to infinity — centered at the origin and placed 
to its right, plus a straight part overlapping the imaginary axis, which avoids the origin, the 
possible purely imaginary zeros ±iyi, yi > 0, and the points ±i/&i, ±i/b 2 by small semicircles 
whose radii tend to zero. However, the contribution from the small semicircle around the origin 
will vanish when we manage to shift the initial <r-domain to the left of its initial position and a 
becomes negative enough to reach — D. Bearing this in mind, we may already neglect this part. 

The asymptotic behaviour of the F function on the upper and lower half-planes motivates 
the factorization 

F(z)=F 1 (z)F 2 (z), 



where 



±1 



F 2 ± (z) 



-{l±ib 1 z){l±ib 2 z)e Tlz for lm(z 

2 (i T *M(iT*M e±2i , forIm(z) > 



(5.6) 



(5.7) 



(l±ibiz)(l±ib 2 z 
On this basis, the original integration is decomposed as follows 

/ dzz-°— In F(z) = [ dzz~ a — In FUz) + / dz z~ c — In F 2 (z) 
Jc dz Jc dz J c dz 



dzz- a — In F 2 (z) 
dz 



dzz- a —\nF+(z)+ [ 
dz ~ Jc 



dz z 



(5.8) 



with C + and C denoting the upper and lower halves of the integration circuit, which have 
Im(z) > and Im(z) < 0, respectively. By virtue of our definitions of F^ 2 , 



In Ff(z) 



0(e ±2iz ) forlm(z)>0. 



(5.9) 



Thanks to these properties, the integrals involving In F^ 1 will vanish on the large circular parts 



as their common radius — say R 

A- stand for Arcs and V* for Vertical,, we have that / 

J A 

and, in consequence, 



tends to infinity, i.e., separating 

d 



dz z 



± 



+ 



A= 



where 



dz 



In F 2 (z) — > 0, as R — ► 00 



1 

1 

" 2vri 



d 



c+ 



L 



dzz- a — In F+(z) + / dzz- a — In FT(z) 
dz 1 



dz 



dz z 



v+ 



dz 



In F+(z) + I dz 
Cf (a) + AC A ' 



dz 



In f; 



(5.10) 
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Here we use the notations 



d°V) 



8111 



p.v. 



dtt~ 



, m=l 



dt 



cos 



In 



(1+M)(1 + M) 
(1 - M)(l - M) 



-2f 



7T<7\ 
V~2~/ 



(5.11) 



Observe that the result denoted by C\\ a ) has been obtained after parametrizing the V ± seg- 
ments in the way: z = e m / 2 t with t from oo to 0, for V + ', and z = e~ tn ^ 2 t with t from to oo, for 
V~ . The A^a(c) term comes from the integrals over semicircles avoiding the purely imaginary 
zeros and poles of F 2 (z). 

The remaining piece in (|5,8|) can be evaluated similarly: 



1 

2ni 



dzz^—ln FAz) 
c dz 



1 7T<7 

— sin — p.v. 

7T 2 F 



dz z 



1 



m=l 



l/6r, 



E 6(bm)bm COS 



m=l 



7T(7 

~~2' 



(5.12) 



where the second summand on the right comes from the poles ±i/b m for b m > 0. Now we see 
that this term cancels the corressponding one in the expression for A£a(<t). Putting together 
all these contributions — eqs. (|5T5|) , Q5.1C| ), ( |5.11| ), ( |5.12| ) — , we arrive at 



A (0), \ 7TCT _ CT 

Q » - cos — 



1 7T<7 

H — sin — p.v. 

7T 2 ^ 



f°° / 2 1 
/ dz z~ a 1 + V — 



(5.13) 



From eqs. ( |5.3| ) and ( 5.13 ) one derives 
e c (/x;s) = 



cotan — ( — ) + p.v 



2 c 7r( D - 1 )/ 2 r(s/2)r(l - a/2) 
1 



2 ^ 

i 



a + 



2 



i * 



d 1 

H In 

dt 



1 



(1 + /?!*)(! + /3 2 t) 



-2ai 



(5.14) 



It can be seen that the contribution to the vacuum energy corresponding to the first term in 
the square brackets under the integral is the vacuum energy density for the Minkowski vacuum 
without boundaries and its regularized value is equal to zero. Hence, we will omit this term in 
the following consideration. The contribution of the complex zeros vanishes at s = — 1 due to 
the factor cotan(-7rs/2). The contribution of the last, logarithmic, term is finite at s = —1 and 
vanishes in the limit a — > oo. It follows from here that the term 



ae 



00 



(/*;*) 



/i s+l 2 -£> 7r (l-£>)/2 



-p.v. 



eft 



D-l-s 



f - 1/& 



(5.15) 



corresponds to the total vacuum energy in the region a m < x < oo for the case of a single plate 
with Robin coefficient (3 m . By taking into account that the regularized value of the integral 
Jq 00 dtt D ~ l is zero we see that E^ s \fj,; —1) coincides with the total vacuum energy derived 
in section [| by summing the corresponding volume and surface energies. The integral in formula 
( [5.151 ) can be evaluated by using the standard formulae (see, for instance, JT9[): 



2 -D 7r (3-D)/2 



;s 



Wm\) 



s+1 



I A, 



l^r(|)r(fi^) sin(D 



S 7T 



1. 

cos(D 



/3 m < 

S)7T, /3 m > 



(5.16) 
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As we see the total vacuum energy has a pole at s 
for the pole structure one derives 



■1. Laurent-expanding where necessary, 



-D-l 



7T 



-D/2 



^r(z>/2 + i) U + i 



1 



+ ln(/u |Ai 



+ -"0 



D 



+ 1 



+ o(s + r 



(5.17) 

where ip{s) = dlnT(s)/ds is the diagamma function. We will assume that the pole term is 
absorbed by the corresponding part in the bare action [jTj. The way in which the pole is removed 
is not unique. The renormalization scheme of ref. |^3| corresponds to the minimal subtraction 



and is equivalent to simply removing the pole term from eq. (5.17). In this way for the finite 
part of the total energy of a single plate we receive 



£ (S) (A 



2hu>|/? m |) + V>(£>/2 + l) -^(-1/2) 
2 D + 2 ir D / 2 p°T(D/2 + 1) 



(5.18) 



This energy depends on the arbitrary scale [/,, as usually happens with this type of regularization. 
The discussion for the role of this scale can be found in |23|| . The corresponding volume and 
surface energies are obtained by using relations ( |3.33|) . Note that the term containing the 
normalization scale is independent on the distance between the plates and does not contribute 
to the vacuum forces. 

-(D 



Now from eq.(5.14), for a 
a\ < x < CL2 one derives 



where 

c-(0) 



1 — s), for the total vacuum energy in the region 



E= Y,E {s) (P m ) + ae 



(o) 

c > 



2 D + 1 ir D / 2 T(D/2 + l)a D + 1 



m=l 



p.v. 



7 

dtt D - In 
dt 



(i + M)(i + M) 

(!-&!*)(! -fct) 



-2t 



(5.19) 



(5.20) 



Further, performing a differentiation in the integral in eq.( |5.20D , one obtains the following de- 
composition 

4 0) = £>+e® , (5-21) 



where e c and e c are defined in ( [4.16 ) and ( 4.32; ) . Now, comparing this with eqs.( 4.25| ), ( |4.16 ), 
we recognize as the volumic part of the integrated energy per unit-volume of a conformally 
coupled field, i.e., the part coming just from the volume between the plates, already calculated 
in the previous section. Since accounts for the "interference" part of the total integrated 



energy per unit- volume, the difference, e 
from the surfaces of the plates. 



(o) (i) 



.(2) 



has to be identified as the contribution 



5.2 Identification in terms surface density 

Next, we may consider the implications of this fact in terms of the densities found in the previous 
section. From relation (4.7), for the integrated Casimir energy per unit volume in the region 
between the plates, one obtains 



(vol) 



1 

a 
1 

2rx 



(12 



dx{0\T 00 \0) 



01 



E 



(2tt) d - 1 



k\ + z 2 /a 2 + 



(4£ — l)z sin2:cos(z + 2a%) 



&2 _|_ z 2 j a 2 1 + \ sin z cos(z + 2a\ ) 



(5.22) 
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where for the purely imaginary zeros the kj_ - integration goes over the region k 2 ^ > yf/a 2 . 
As we see, this result differs from the total Casimir energy per unit volume ( |5.1| ). We have 
argued that the reason for this difference should be the existence of an additional surface energy 
contribution to the volume energy (5.22), located on the boundaries 1,2. The 

corresponding energy density is defined by relation (|3.29| ), where now 



5(x; dM) = 5(x - a x + 0) - 6(x - a 2 - 0). 



(5.23) 



The corresponding v.e.v. can be evaluated by the standard method explained in the previous 
section. This leads to the formula 



(0 



p(surf) 
00 



0) 



4£- 1 
2 

4g- 1 



S(x;dM) (d x {0\<p(xMx')\0)) 
5(x;dM) ' ' 



E 

z=\„,iyi 



sin [2z^ + 2ai] 



(27T) 



D-l 



_|_ z 2 j a 2 1 + \ sinzcos(,z + 2a\) ' 

(5.24) 

which provides the energy density on the plates themselves. Then, for the integrated surface 
energy per unit area one obtains 

1 



(surf) 



dx (0 T o f f) 



4£- 



1 



2a 3 



E 

z=\ n ,iy l 



0) 

d D - 1 k l 



sinzcos(2; + 2ot\) 



(5.25) 



(27T) 



D-l 



^2 _|_ z 2^ a 2 1 + j sinzcos(z + 2ai] 



Adding up ( 5.22| ) and ( |5, 25 ), one re-obtains the standard result (|5.1|). Thus, we have just checked 
that e, 

may be written as follows 



£ (voi) _j_ £ ( sur fT^ After the standard integration over transverse momentum eq. ( 5.25 ) 



.(surf) 



4£-lDT(-D/2) 

a D+l 2 D+l n D/2 



E : 

z=\„,±iyi 



D 



z=X n ,±iyi 



+ - sinzcos^ + 2a\) 



(5.26) 



The regularized value for the first sum in the square brackets have been found in the previous 
subsection (see ( |5.11 ) and ( |5.13| )). The second sum might be evaluated using the Abel-Plana 
summation formula. Applying formula (B^) to this sum and omitting the divergent contribution 
from the integral term (this corresponds to the calculation of (5.24) taking in this formula 

(<p(x)<p(x'))suB = (Q\v(x)<p(x')\0) - (0m\<p(x)<P&)\0m) instead of (0 \ip(x)<p(x')\ 0)) for the 
surface energy per unit area one obtains 



^(surf) 



as 



(surf) 



£150 

m=l 



s)(surf ) 



(Pr, 



aL>(4£ - 1)4 2) 



(5.27) 



with e ( c > defined in flOp . Now using f[2l|) and (gjgg) we can find the volume part of the 
vacuum energy as 

2 

E {vo\) = E (tot) _ ^(surf) = £ E ( B )(vol) (/?m) + a Ul) + 4jD( £ _ ^) £ (2)\ ^ (5 _ 2g) 



m=l 



which coincides with the previous result (4.31) obtained by integrating energy density. Hence, 
we have shown that the local and global approaches lead to the same expression for the volume 
energy in the case of the scalar field with general conformal coupling £. Note that, as it follows 
from fl5.27p , the quantity Sc is the additional (to the single plate) surface energy per unit 
volume in the case of the conformally coupled scalar field. 
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5.3 Numerical examples 



Formulas (5.20) and ( 5.2 1| ) with ( |4.16j ), ( |4.32j ) are suitable for numerical evaluation (actually, we 



have used an alternative form of found by partial integration of the integral in eq.( |5.20| )). 
Making use of these expressions, we may, e.g., keep the value of b\ fixed and study the varia- 
tion with &2 of the total — — , volumic — Ec — , and superficial — e^f 1 — integrated Casimir 
energies per unit-volume in the conformal case. An example for D = 3 is given in Fig. ||, 
where we have set 6i = while 62 changes. The curve in solid line depicts a^ef^ and the 
one in dotted line a 4 £^\ being the surface contribution the difference between them, 

which has been plotted in dashed line. Note that, at 62 = 0, one recovers the known result 
a 4 e = — jjr^ T(2) Cr(4) = — 1^40 — ~~ 0.0069 for the total and volumic parts, while the super- 
ficial contribution is zero. The volume part is higher than the total result, meaning that the 
surface contribution is always negative. In fact, the magnitude of the latter tends to zero when 
62 — ► —00, as had to be expected, because this limit corresponds to Neumann boundary condi- 
tions. Thus, between 62 = and that asymptotic regime it must have at least one minimum, 
and we actually observe one at 62 — —0.70. To be remarked is the zero of 4 0) at 6 2 ^ -0.81. 
Below this value, the total Casimir effect is repulsive while, for larger values of 62 until 62 = 0, 
is attractive. 

It is also possible to consider the Casimir energy as a function of b±, 62 simultaneously. The 
plot shown in Fig. illustrates the changing form of the total energy e c = ef^ on a given region 
of the (61, 62)-plane (for D = 3, too). An analogous description of e c is provided by Fig. |5[ 




Figure 3: Total integrated Casimir energy per unit-volume multiplied by a 4 , i.e. a 4 e c , for 
D = 3, bi = 0, and — 5 < 62 < 0. Separately plotted are the total value, in solid line, the 
volumic contribution, in dotted line, and the superficial part, in dashed line. Note the minimum 
of this surface contribution at 62 — —0.70, and the zeros of the total value, at 62 — —0.81, and 
of the volume part, at 62 — —0.58. 



6 Ending comments 

In the present work we have dealt with a calculation of the Casimir energy when one sets 
Robin boundary conditions on one single plate or a pair of parallel plates. Its evaluation has 
been based on a variant of the generalized Abel-Plana summation formula in ref.[|l|], adapted 
to these situations, and derived in the appendix B. This method turns out to be adequate for 
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Figure 4: Left: Total integrated Casimir energy per unit-volume times a 4 , i.e. a 4 s c , correspond- 
ing to D = 3, —2 < b\ < 0, —2 < 62 < (in this region, there are no imaginary eigenfrequencies 
and the p. v. prescription for the integrals is unnecessary). Right: contour representation of the 
same plot. 




20 



finding vacuum expectation values of the energy-momentum tensor, i.e., local densities. From a 
slightly different viewpoint, zeta function regularization has been applied to the summation of 
eigenfrequencies, which directly gives the integrated energy per unit-volume. 

When just one plate is considered, the only present parameter is the relative coefficient 
between the non-derivative and derivative terms in the boundary condition The local 

density is given by formula ( p. 19 ), which vanishes for the conformal value of the curvature 
coupling. Otherwise, this formula depicts the local dependence of this density (exemplified by 
Fig.|l] for D = 1 and D = 3), which is singular on the plate itself. Note that the requirement 
of conformal invariance has the power of suppressing the presence of divergent parts, just as 
happened — for a different system — in ref . || . 

If there are two parallel plates, the relevant parameters are three: the (rescaled) relative 
coefficients between the non-derivative and derivative parts at each boundary (pi and 62); an d 
the separation length between them (a). Then, the additional (to a single plate) total integrated 
Casimir energy per unit-volume is given by formula ( |5.20D , and its decomposition into purely- 
volume and purely-surface parts by eqs. ( |4.31| ) and ( |5.27 ), in terms of the quantities and 
e^ 2 ) defined by eqs. ( 4.16| ) and ( 4.32 ). If the coupling is conformal (£ = £ c ), e^ 1 ^ and 
themselves coincide with the volume and surface contributions, respectively, and, in any case, 



the decomposition ( 5.21 ) holds. The surface contribution, coming from the plates themselves, 
would be absent for Dirichlet or Neumann boundary conditions. 

To be remarked is the fact that, at least in some situations free of imaginary eigenfrequencies, 
there are parameter choices which give a vanishing Casimir energy. As illustrated by Fig. ||, one 
may vary the value of the 62 parameter so as to reverse the sign of the effect. At the same time, 
we have seen that there is another Rvalue for which the surface contribution has a minimum. 
Examples of simultaneous variations of bi and 62 are shown in Figs. |4] and ||. 

An interesting feature of the Casimir effect with Robin boundary conditions is that there is 
a region in the space of parameters defining the boundary conditions in which the vacuum forces 
are repulsive for small distances and attractive for large distances. This leads to the possibility 
for the stabilization of the plates separation by using the Casimir effect. 
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A Appendix: Complex zeros 



First of all we will show that the real and possible purely imaginary zeros (see below) of F(z) are 



simple. To see this, we note that on the class of solutions to (4.4), the corresponding derivative 
can be presented in the form 



F'(z) 

Using the integral relation 



1 



1 H — sin z cos(z + 2a,\ 

z 



[(l — 6i&2^ 2 ) cos z + (b% + b\)z sin , 



1 f 1 

lH — sin z cos(z + 2a\) = 2 / dx cos 2 (zx + a\) 
z Jo 



(A.l) 



(A.2) 



we conclude from here that F'(z) 7^ if z, z 7^ is a zero of F(z), and hence these zeros are 
simple. 
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Purely imaginary zeros of F(z) may exist. This sort of solution has to do with the presence 
of imaginary parts in the eigenfrequencies. They can be detected as the real zeros of the de- 
nominator in the last integral of eq.(B.4). It is convenient to write the corresponding equation 
in the form 

(bx + b 2 )t , , 

tanht = ; , , V A.3 

After studying the nature of this equation in terms of b\ and b 2 , one finds out that: 
1) Equation ( |A.3| ) has no positive real zeros for 



{61 + 62 > 1, 6i6 2 < 0} U {b 1>2 < 0} . (A.4) 

2) Equation ( |A.3| ) has a single positive real zero for 

{0 < 61 + 6 2 < 1, hb 2 < 0} U {61 + 6 2 > 1, 6i,a > 0} U {b x + b 2 < 0, b x b 2 < 0} . (A.5) 

3) Equation ( |A.3| ) has two positive real zeros for 

61 + b 2 < 1, &i )2 > 0. 

The parameter values in Fig. ||, namely b\ = and —5 < b 2 < 0, fall into case 1), when 
there are no real positive zeros. When F(z) has complex zeros (situations 2 or 3), their extra 
contribution to the mode sum is given by formula (|B.8|) in appendix B. 



B Appendix: Summation formula 

The vacuum expectation values for the physical quantities in the region between plates will 
contain the sums over zeros of the function F(z) defined by (|4,4|). To obtain the summation 
formula over these zeros we will use the generalized Abel-Plana formula (GAPF) ffl. In this 
formula as a function g(z) let us choose 

f(z) 

g(z) = -i [(I-6162Z 2 ) cos z + (b 2 + h)z sin z] (B.l) 

b \z) 



For the sum and difference in the GAPF one has 



g(z) ± /(*) = i(b lZ ± i)(62^ ± ^ ±IZ Tt\- (B.2) 

F{z) 

Let us denote by A n , n = 1, 2, . . . the zeros of the function F{z) in the right half-plane, arranged 
in ascending order, and by ±iyi, yi > the possible purely imaginary zeros of this function. It 
can be easily seen that 

Res z=An g(z) = — — ~ %f{z) (B.3) 
1 + - sin z cos(z + la\) 



(as it follows from ( A. 2 ) the denominator on the right of this formula is always positive). First, 
we will consider the case of function f(z) analytic for Rez > 0. Now substituting in GAPF 
(formulas (2.10)-(2.11) in [Q]) (BJJ, ( |B.2 ), taking the limit a — > (here a is the parameter on 



the right of GAPF and the poles Hyi are excluded by small semicircles with radius p on the 
right half plane, p — > 0), and using (|B.3[) one obtains the following summation formula 

v ' ^ = -* m + rnz)dz 

^ l + isinzcos(z + 2ai) 2 1 - b 2 - bx 



=\ n ,±iy. 







/ [f(it)-f(-it)]dt 
(6 1 t+l)(6 2 *+l) 
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where the prime on the summation sign means that the contribution of terms corresponding 
to the purely imaginary zeros have to be halved. This contribution comes from the integrals 
taken around semicircles enclosing these zeros. Note that the denominator on the left can be 
also written in the form 

1 + - smzcoBfr + 2a!) = 1 - (1 + + ^ , z = X n ,± m . (B.5) 

In flB.4| ) we have assumed that b\ + b 2 ^ 1. In the case b\ + b 2 = 1 to ensure the convergence at 
origin in the second integral on the right of formula ( |B.4[) we need to have f{z) ~ faz a , a > 2, 
z — ► 0. Now the first summand on the right of this formula should be replaced by 

nfo6a2 (B.6) 



2(b( -h + 1/3) 



Formula (R4) is valid for functions f(z) satisfyng the condition 

\f(z)\ < e{x)e clyl , z = x + iy, \z\ -> oo, (B.7) 

where c < 2, e(x) — » oo for x — * oo, and having no poles on the imaginary axis. However, as 
follows from (|4.10|) , ( |4. 1 1|) , for a scalar field with C 7^ Cc the corresponding function has the form 
f(z) = f„ l {z)/{z 2 b 2 n + 1), m = 1,2. In this case the subintegrand on the right of GAPF has 
purely imaginary poles ±i/b m for b m > 0, m = 1, 2. In analogy to the purely imaginary zeros of 
F(z), these poles have to be excluded from the integral over the imaginary axis by semicircles 
on the right-half plane. The integrals over these semicircles will give additional contributions 



7T 



2 \ b m \ 



e(b m ) [f m {e m/2 / \b m \) + f m {e-™ /2 / \b m \)] . (B.8) 

to the right-hand side of flB,4| ). In the case b\ = —b 2 and for functions f{z) having no poles on 
the imaginary axis from (|B.4D one obtains the Abel-Plana formula in the usual form. 
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